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ABSTRACT 

The gravitational lens system CLASS B2 108+2 13 has two radio-loud lensed images sepa- 
rated by 4.56 arcsec. The relatively large image separation implies that the lensing is caused 
by a group of galaxies. In this paper, new optical imaging and spectroscopic data for the lens- 
ing galaxies of B2 108+2 13 and the surrounding field galaxies are presented. These data are 
used to investigate the mass and composition of the lensing structure. The redshift and stel- 
lar velocity dispersion of the main lensing galaxy (Gl) are found to be z = 0.3648 ± 0.0002 
and (Ty = 325 ± 25 km s~^, respectively. The optical spectrum of the lensed quasar shows no 
obvious emission or absorption features and is consistent with a BL Lac type radio source. 
However, the tentative detection of the G-band and Mg-b absorption lines, and a break in the 
spectrum of the host galaxy of the lensed quasar gives a likely source redshift of z = 0.67. 
Spectroscopy of the field around B2 108+2 13 finds 51 galaxies at a similar redshift to Gl, 
thus confirming that there is a much larger structure at z ~ 0.365 associated with this system. 
The width of the group velocity distribution is 694 ± 93 km s~^, but is non-Gaussian, imply- 
ing that the structure is not yet viralized. The main lensing galaxy is also the brightest group 
member and has a surface brightness profile consistent with a typical cD galaxy. A lensing 
and dynamics analysis of the mass distribution, which also includes the newly found group 
members, finds that the logarithmic slope of the mass density profile is on average isothermal 
inside the Einstein radius, but steeper at the location of the Einstein radius. This apparent 
change in slope can be accounted for if an external convergence gradient, representing the 
underlying parent halo of the galaxy group, is included in the mass model. 

Key words: gravitational lensing - quasars: individual: CLASS B2108+213 - cosmology: 
observations. 



1 INTRODUCTION 

In the standard model of galaxy formation, structure arises hierar- 
chically through the mergers and interactions of smaller mass sys- 
tems, leading to increasingly more massive structures. Cold dark 
matter (CDM) simulations of these merging events predict uni- 
versal radial matter density profiles with inner slopes of 7 ^ 1- 

which are indepe ndent of the 



1.5 [where p cDM(r) oc 

halo niass [e.g.lNavarro. Frenk & W hite 199 ^ (NFW);lMoore et all 
ll998l : lDiemand. Kuhlen & Madaull2007i : [Navarro et aL"201(y]. Fur- 
thermore, the compact cores of merging systems are expected 
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to survive the cannibalism pr ocess, forming CDM substructure 
within the new paren t halo (e.g. lMoore et al.l 19991 : iGaoet al.l2004l : 
iDiemand et"aDl2007h . 

Strong gravitational lensing, where multiple images of the 
same background object are formed by an intervening potential, 
has provided a powerful observational test of these theoretical pre- 
dictions for mass structures out to 2; ^ 1, particularly when com- 
bined with stellar kinematic data (e.g . iTreu & KoopmansI l2004l : 
ISand et al.ll2004l : lKoopmans et al.ll2006l) . 

In the case of galaxy-scale lensing (Miens =10^^- 
10^^ Mo), detailed mass modeling of the positions and flux 
densities of the multiply imaged source has shown the to- 
tal (i.e. baryons plus dark matter; ptot) mass density pro- 
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files of lens galaxies to b e close to isothermal, that is, 
7^ = d\o% ptot/dlo^r ^ 2 fcohn et aDl200ll: Koop mans et alj 
20061: iRusin. Kochanek & K eeton 200 31: IXreu & K oopmans 200i; 



Wucknitz. Big gs & Bro wnd 12004 : iKoopmans et al.! i2009) with 



very little scatter (^5 per cent on 7^). The steeper slope than that 
predicted by dark matter only simulations is generally attributed to 
baryon cooling and collapse, which steepens the inner density pro- 
file for^lens g alaxies (|Blumenthal et al. 1984; Gnedin et al. 2004; 
lKeeton"2001b: Kochanek & Whit e.2001i : [Jiang & Kochanek 2007^). 
The origin of the small scatter in the total mass density pro file slope 
is not completely clear (see discussion in lKoopmans et al.ll200d) . 

Decomposing the baryonic and dark matter components is par- 
ticularly hard for galaxy-scale lenses, given the dominance of bary- 
onic mass at small radii. By combining lensing with stellar kine- 
matic tracers, and assuming a spatially uniform mass-to-light ratio 
for the stellar distribution, an upper limi t of 7 ^ 1.3 for the dark 
matter component of galaxy- scale haloes jTreu & Koopmansll2004l) 
can be obtained, which is certainly within the allowable range pre- 
dicted from simulations. However, studies of the inner slopes of 
dark matter haloes using non-lensing techniques, for example, from 
rotation curve measurements of nearby dwarf and spiral galaxies, 
find evidence for a flattened core (0 < 7 < 1) out as far as the 
optical rad ius (e.g. Salucci & Burkert I2OOOI : IS waters etaDl2003l : 
iGentile et al.ll2004l : lDonato et al.ll2009l) . which are typically incon- 
sistent with the CDM models. Although, it should be noted that 
due to systematic effects, such as non-circular motions or poorly 
defined disk i nclinations, the CDM models cannot be discounted 
jSwaters et al . 2003). 

For galaxy clusters (Miens > 10^^ M©), the situation is sig- 
nificantly different. On the one hand, the stellar component is less 
important than for galaxies on the scales that are relevant to mea- 
sure the inner slope of the dark matter, and several mass probes 
are available to cover the broad range of radii needed for this mea- 
surement, including strong lensing, stellar kinematics. X-rays and 
weak lensing. On the other hand, many clusters are not dynamically 
relaxed systems, often with significant triaxiality and substructure. 
Therefore the interpretation of the results can be more complicated 
than in the case of galaxies and can be affected by potential system- 
atics, such as those arising from projection effects and departures 
from hydrostatic equilibrium. 

A relatively large scatter of inner de nsity profi le slopes 
has been reporte d for galaxy clusters ( Smith et al.l | 200l |: 
Sand. Treu & ElHsl I2OO2I; lEttori. De Grandi & Molendil l2002l: 
Lewis. Buote & Sto cke 20031: ISandetal.1 |2004 l2005l: bavazSl 
2OO5I: ISaha. Read & Wilhalm l2006l : ISharon et all l2005l : 
Bradac e t al. 2008) - ranging from flat inner cores to slopes 



consistent with CDM predictions. However, it should be noted 
that some of these authors refer to the slope of the dark matter 
component (7) or to the slope of the total mass density profile (7O 
as if they were the same quantity. This is not justified. In fact, 
although stars are only a small fraction of the total mass, they are 
extremely spatially co ncentrated and the refore they affect the slope 
in a significant way (' Sand et all 120021) . Overall, although more 
work remains to be done at cluster scales, it seems clear that in at 
least some cases the inner slope of the dark matter halo is flatter 
than predicted by CDM simulations and the NFW profile is not a 
good fit to the data dSand et al.ll2008h . 

Galaxy groups (Miens = 10^^-10^^ M©), containing a few 
to a few tens of galaxies, are the mos t common envir onment for 
galaxies to reside at low redshifts fe.g. lEke et al.ll2004b . Although 
the total mass of a galaxy group can be estimated from the luminos- 
ity distribution or the stellar mass of the individual group members 



(e.g.lYang e t al.'2007*) a s well as from the group velocity dispersion 
( Zabludoff & Mulchaevll 19981) . the distribution of the matter within 
groups is not well known. 

Galaxy groups have been relatively unexplored with strong 
gravitational lensing due to the paucity of known direct im- 
age s plitting by grou p sized haloes (see iLimousin et al ] l2009al 
and tKubo et al.l l2009l for recent developments) and weak lens- 
ing measurements can require the stacking of possibly hundreds 
to th ousands of groups (e.g. Hoe kstra et al.l EoOll : IParker et aP 
I2OO5I : .Mandelb aum et al.l [2OO6), and even then, only probe the 
outer part of the dark matter halo. There are several examples 
where strong gravitational lensing galaxies have been found to 
be associated with a group of galaxies and/or there are groups 
of galaxies found along the line-of- sig ht to the lensed o bject 
(Tonry 1998: T onrv & Koc hanek 19 99. 200QI: Ipassnacht & LubinI 
2002l?lFaure eTal . 2004; M organ et al . 2005; F assnacht et al.ll2006l : 
Momcheva et al.ll2006 : .Wimams et al.>.2006 : .Auger et al.ll2007h . In 



these cases, the lensed image separations are still in the galaxy- 
scale regime of ^0.5-1.5 arcsec because the lensing potential is 
dominated by a single galaxy, with the contributi on of the group to 
the overall convergence being only ^5 per cent jMomcheva et al.l 
I2OO6I : lAugeretal.ll2007l . l2008h . However, it should be noted that 
environmental contributions at this level can still have a large 
impac t on the results of lens modelling ( Keeton & Zabludora 
I2OO4I) . For lens systems with image separations larger than 3 arc- 
sec the contribution of the environment is thought to be much 
larger than in the galaxy- scale regime (with an external conver- 
gence > 10 per cent), p ossibly enhancing the lensing probabilities 
JOguri. Keeton & DalafcoOSh . 

How the mass is distributed within a group- sized halo, and 
what effect baryons have on the overall density profile is not en- 
tirely clear. However, as the image separation increases there is ex- 
pected to be a transition from galaxy scale haloes, which are dom- 
inated by stars in the centre and have a relatively steep total mass 
density profile (i.e. 7^ ~ 2), to group-scale haloes, characterized by 
a flatter inner density profile of 7^ ^ 1-1.5 (|Oguri 2006). Another 
question to be addressed is whether the dense environment affects 
the dark matter sub-haloes of the satellite galaxies within the group. 
It is expected that due to tidal interactions, the dark matter of the 
sub-haloes will be strip ped, res ulting in a change in the sub-halo 
mass-to-light ratio (e.g. iGao et al. 2004) and a truncation which 
would lead to a mo re concentrated sub-halo mass distribution (e.g. 
iBullock et all200lh . 

In this paper, we present for the first time an analysis of the to- 
tal (baryonic and dark) inner density profile of a moderate redshift 
galaxy group using strong gravitational lensing and stellar kine- 
matics. We show that the gravitational lens system B2 108-1-2 13 is 
part of a massive group of galaxies and that the multiple imaging 
is found around the large central galaxy, as opposed to the satellite 
members of the group. Thus, we are probing the inner part of the 
group-scale halo. In Section |2l we present an introduction to the 
lens system and summarize the previous observations. In Section 
[3] new ground based imaging with the W. M. Keck Telescope is 
combined with archival HST data to select potential members of 
the galaxy group associated with B2108-I-213. Our spectroscopy of 
group candidates is presented in Section |4] From these observa- 
tions we found 5 1 galaxies at the same redshift as the main lensing 
galaxy. The properties of the group and the constituent group mem- 
bers are discussed in Section |5] An analysis of the stellar velocity 
dispersion and surface brightness profile of the main lensing galaxy 
is presented in Section[6] In Section|7]we present a new mass model 
for B2 108-1-2 13 which incorporates the group galaxies and the stel- 



© 2008 RAS, MNRAS 000,[T}{T8] 



Gravitational lensing by a moderate redshift galaxy group 3 



lar kinematics. We discuss our results in Section [8] and present our 

conclusions in Section |9] 

For all calculations we adopt an Qm = 0.3, 
=0.7 spatially-flat Universe, with a Hubble constant of 

Hq — 70 kms~^ Mpc~^. At redshift 0.365, 1 arcsec corresponds 

to an angular size of 5.075 kpc. 



2 CLASS B2108+213 

The gravitational lens system CLASS B2108+213 (iMcKean et alJ 
[2005.) . which was discove red in the final pha s e of the Cosmic Len s 
All-Sky Survey rCLASSi lM^ers et al.ll2003l : lBrowne et al.ll2003h . 
consists of at least two lensed images of a radio-loud quasar that 
are separated by 4.56 arcsec. The wide image separation between 
the two lensed images implies lensing by a massive halo, due to 
a group-sized mass distribution. Therefore, B2 108-1-2 13 provides 
an excellent opportunity to probe a mass regime where very lit- 
tle is known. The flux-ratio of the two lensed images (A and B; 
where image A has the higher flux-density) is ^2:1. High resolu- 
tion radio imaging at 1.7 and 5 GHz with Very Long Baseline In- 
terferometry (VLB I) found a compact core and extended jet emis- 
sion in both lensed imag es, which is consistent with gravitational 
lensing ( More et al. 2008). From optical and infrared imaging with 
the Hubble Space Telescope (HST) sl massive early-type elliptical 
galaxy (Gl) was found between the two lensed images, which is 
thought to dominate the lensing potential. However, there is an ad- 
ditional satellite galaxy (G2) within the Einstein radius and a scat- 
ter of nearby galaxies of similar colour to Gl within ^15 arcsec 
of the system. In addition, the infrared emission from lensed image 
A shows a large gravitational arc which is due to th e lensed emis- 
sion f rom the host galaxy of the background quasar (iMcKean et al.l 
l2005h . Observations with Chandra revealed luminous X-ray emis- 
sion, which is presumably from the hot group/cluster gas, that is 
elongat ed by ^ 1 Mpc in an east-west direction across the lens 
system (iFassnacht et al.|[2008h . There is also a third radio compo- 
nent (C) which is coincident with the optical emission from Gl to 
within the astrometric uncertainties. However, component C is not 
thought to be a third lensed image because, first, further radio imag- 
ing with VLB I and the Multi-Element Radio Linked Interferome- 
try Network (MERLIN) found compact and extended radio emis- 
sion consistent with an Active Galactic Nucleus (AGN) embedded 
within the main lensing galaxy Gl, and second, the flux-density 
of component C is more than 10^ times too high when compared 
to what is expect ed from modelling the lensing mass distribution 
(iMore et al.ll2008h . 

A simple lens model, which includes both Gl and G2 with 
isothermal density profiles and an external shear, can explain the 
positions of the lensed images A and B, bu t fails to reproduc e 
the observed flux-ratio (iMcKean et aP l2005l : [More et al.l l2008h . 
This could be explained if there is a perturbation in the lensing 
mass model, for example, due to mass substructure in the parent 
dark matter halo (e.g. Bradac et aL 20021: IPalal & Kochanekll2002l : 
iMetcalf & Zhaoll200ll\^getti et aLll201Ql) or from the contribution 
of the nearby group galaxies. The best fitting model for the system 
that satisfies the positional and flux constraints requires the com- 
bined luminous and dark matter density profile of Gl to be steeper 
than isothermal ( 7^ = 2.451q; at the 68 per cent confidence level; 
iMore et al.ll2008h . 



Table 1. Observing logs for the photometric and spectroscopic (long- and 
multi-slit) observations of the B2 108+2 13 field. These observations were 
carried out with the W. M. Keck-I Telescope. Slit position angles (PA) are 
given east of north. 



Filter or slit/ Date Exposure time (s) PA 

mask number LRIS-B LRIS-R (°) 

U 
B 
V 
R 
I 

Slit 1 
Slit 2 
Slit 3 
Slit 4 
Slit 5 
Mask 1 
Mask 2 
Mask 4 
Masks 
Mask 6 
Mask? 



3 PHOTOMETRY 

The first step of our investigation into the local environment and 
lensing galaxies of B2 108-1-2 13 was to obtain deep imaging of the 
field using the W. M. Keck Telescope. These observations were 
designed to identify galaxies which could potentially be part of a 
group or cluster associated with the main lensing galaxy, Gl. 

3.1 W. M. Keck Telescope imaging 

We carried out extensive optical imaging of the B2 108-1-2 13 field 
with the W. M. Keck-I Telescope on 2003 October 24 and 25. The 
conditions on both nights were not photometric due to faint cirus 
cloud cover. Therefore, photometric redshifts of the field galaxies 
could not be calculated. However, the relative colours could still be 
used to select group candidates for follow-up spectroscopy. The full 
width at half maximum (FWHM) seeing was ^0.7 arcsec through- 
out the observing run. Our imaging data were taken with the Low 
Resolution Imaging Spectrograph (LRIS; lOke et al.llT99? ). which 
allows simultaneous observations through blue (-B) and red (-R) 
optimized cameras. A 5600 A dichroic was used to split the light 
between LRIS-B and -R for our U, B, R and /-band filter photome- 
try. For our V-band imaging, we employed a 4600 A dichroic. With 
these filters in place, the usable field of view through LRIS-B and 
-R is ^ 5.5 x7.5 arcmin^. The integration times through each filter 
were typically n x 300 s, with a 10 arcsec dither between each ex- 
posure to aid the removal of bad pixels and cosmic rays during the 
reduction stage. A summary of our imaging observations is given 
in Table[T] 

The data were reduced within lRAl(3 using standard proce- 
dures. The world coordinate system was added to each image by 
fitting to the known positions of United States Naval Observa- 
tory (USNO) stars. This resulted in an astrometric uncertainty of 

^ IRAF (image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatories, which are operated by AURA, 
Inc., under cooperative agreement with the National Science Foundation. 



2003 Oct 24/25 


3400 




334 


2003 Oct 24/25 


2400 




334 


2003 Oct 24/25 




600 


334 


2003 Oct 24/25 




2000 


334 


2003 Oct 24/25 




2400 


334 


2003 Oct 24 


5400 


5400 


142 


2003 Oct 25 


2500 


2500 


60 


2003 Oct 25 


2400 


2400 


25 


2003 Oct 25 


1800 


1800 


75 


2007 Oct 10 




2400 


50 


2004 Jun 15 


3600 


3600 


334 


2004 Aug 12 


5400 


5400 


334 


2004 Aug 12 


3600 


3600 


334 


2005 Aug 1 


5400 


5400 


345 


2005 Aug 1 


5400 


5400 


345 


2005 Aug 1 


5400 


5400 


345 



© 2008 RAS, MNRAS 000,[T}tT8] 



4 /. P. McKean et al 



^-^0.3 arcs ec. The photometric ca talogue was generated by SEX- 
TRACTOR teertin & Arnoutj[T996 ). 

No zero-point magnitude caHbration or atmospheric extinc- 
tion correction was appUed to the data because the observing con- 
ditions were not photometric. However, we estimated the apparent 
magnitudes of the objects in the field by bootstrapping an approx- 
imate calibration to the B and R band imaging using unsaturated 
stars which were detected in the USNO Bl catalogue. We used 
archival //5r imaging at F555W and F814W (see Section [T2l ) to 
apply a photometric calibration to the V and / band datasets. 

3.2 Hubble Space Telescope imaging 

HST imaging of the B2108-h213 field with the Advanced Cam- 
era for Surveys (ACS) and the Near Infrare d Camera/Multi - Objec t 
Spectrograph (NICMOS) was presented bv lMcKean etaP ( l2005h . 
To summarise, B2 108-1-2 13 was observed with the ACS through 
the F555W and F814W filters. Both observations were carried out 
in the wide field channel mode, which provided a field of view of 
202 X 202 arcsec^. An infrared image at F160W was aquired with 
NICMOS using the NIC2 camera. However, infrared information 
is only available for 12 sources nearby to B2 108-1-2 13 because the 
NIC2 field of view is 19 x 19 arcsec^. The photometric analysis of 
the ACS and NICMOS data was carried out with SEXTRACTOR. 



4 SPECTROSCOPY 

Spectroscopic observations of B2 108-1-2 13 and the candidate group 
galaxies were carried out with LRIS on the W. M. Keck-I Telescope 
in both long- and multi-slit modes. 

4.1 Long slit spectroscopy 

The long slit data were taken with LRIS during two observing 
sessions. The first set of spectroscopic observations were carried 
out at the same time as the LRIS photometry presented in Section 
13.11 (2003 October 24 and 25). In total, we employed four 1 arc- 
sec wide and 120 arcsec long slits across the field (see Fig. [T] for 
the slit positions). Slit 1 was positioned and orientated to include 
the B2 108-1-2 13 lensing galaxies and lensed images. The three re- 
maining long slits (Slits 2 to 4) targeted those galaxies with similar 
colours to Gl (based on an initial analysis of the LRIS imaging 
data done at the telescope) and in close proximity to B2108-I-213. 
The second set of spectroscopic observations were carried out du- 
ing 2007 October 10. Here, a single 1 arcsec- wide long slit (slit 
5; see Fig.[T]) was placed across the gravitational arc of lensed im- 
age A. The aim of this observation was to measure the redshift of 
the background source. A summary of the total exposure times and 
position angles used for each slit is presented in Table [T] 

For the observations taken on 2003 October 24 and 25, the 
light was split between LRIS-B and -R using a 5600 A dichroic. A 
600 lines mm~^ grism blazed to 4000 A was used through LRIS- 
B. This produced a dispersion of 0.61 A pixel ~^ and a wavelength 
coverage of ^3400 to 5600 A. The light was dispersed through 
LRIS-R using a 600 grooves mm~^ grating centred at 7050 A. This 
resulted in an LRIS-R dispersion of 1.26 A pixel" ^ and a wave- 
length range of ^-^5660 to 8250 A. The wavelength calibration was 
established using Hg, Cd, Ne and Ar arc exposures. An observation 
of the spectrophotometric standard star BD-h28°4211 (IOkelll99Ql) 
was taken on each night to correct for the response of the instru- 
ment. For the observations taken on 2007 October 10, the 5600 A 



dichroic was used with the 400 grooves mm~^ grating centred at 
8500 A through LRIS-R. This gave a dispersion of 1.86 A pixel"^ 
and a wavelength range of ^5700 to 9200 A. Note that no standard 
star observation was carried out so the resulting spectrum of lensed 
image A was not flux-calibrated. The one-dimensional spectrum of 
each detected obje ct was reduced within IRAF following the proce- 
dures described in lMcKean et al.l (|2004|) . 

Our long slit observations yielded spectra for 16 galaxies and 
2 stars. Eight galaxies (including Gl) were found to be at a red- 
shift of ^0.365; thus confirming that the main lensing galaxy of 
B2 108-1-2 13 is part of a larger structure. A detailed discussion of 
the results from our long slit spectra is given in Section[5] 

4.2 Multiple slit spectroscopy 

Our long slit observations only sampled a small number of galaxies 
within --0.3 Mpc of B2108-h213. However, the full extent of the 
group or cluster associated with the B2 108-1-2 13 lensing potential 
may be much larger. Therefore, we undertook multiple slit spec- 
troscopy of the field using LRIS on 2004 June 15, 2004 August 12 
and 2005 August 1. 

The multiple slit capability of LRIS allows the spectra of 
^25-30 objects to be obtained simultaneously within a 8 x 6 
arcmin^ field through the use of slit masks. Candidate group galax- 
ies were selected primarily by colour. In Figs.|2]and[3]we show the 
V—R and R—I colour-magnitude diagrams with respect to R and / 
for the B2108-h213 field. We stress that these LRIS imaging data 
were taken on non-photometric conditions and as such the absolute 
colours provide limited information about a galaxy's redshift. How- 
ever, it is still possible to use the relative colours to identify galaxies 
at similar redshifts. The 5 closest group members to the lensed im- 
ages discovered from our long- slit observations have been labelled 
with red diamonds in Figs. |2] and (3] there is clear evidence of a 
red sequence in both color-magnitude diagrams. On Figs. |2] and [3] 
we have also defined the V—R and R—I colour boundaries which 
were used to select candidate group galaxies. These colour bound- 
aries were set to include all of the group galaxies found from the 
long- slit data, with the exception of the bluest and reddest galaxies. 
Targets were chosen for the slit-masks using four priority classes; 
that is, galaxies meeting i) both the V—R and /colour ranges, ii) 
only the V—R colour range, iii) only the R—I colour range, and iv) 
neither the V—R and R—I colour ranges. The final class was used 
to ensure as many spectra as possible were taken on each slit mask. 

The first set of multiple- slit data were taken during poor 
observing conditions on 2004 June 15. The 400 lines mm~^ 
grism blazed to 3400 A (1.09 A pixel" ^ dispersion) and the 400 
lines mm~^ grating (1.86 A pixel" ^ dispersion) centred to 6800 A 
was used through LRIS-B and -R, respectively. The light was split 
with the 4600 A dichroic. Spectroscopy with a single slit mask 
(mask 1) for 1 h targeted 32 galaxies in the B2 108-1-2 13 field. The 
second set of observations were carried out on 2004 August 12, 
again during poor observing conditions. The 600 lines mm~^ grism 
blazed to 4000 A (0.63 A pixel" ^ dispersion) was used with LRIS- 
B, and the 600 lines mm~^ grating (1.28 A pixel" ^ dispersion) cen- 
tred to 7000 A was used through LRIS-R. Two slit masks (mask 2 
and mask 4) targeted 60 galaxies. The total observing time for each 
mask was 1.5 and 1 h, respectively. The third set of multiple-slit ob- 
servations were carried out during clear conditions on 2005 August 
1 . Here, three slit-masks were employed, which targeted 77 galax- 
ies in the B2 108-1-2 13 field. We also placed a single slit over both 
G2 and lensed image A in each of the three slit masks (total expo- 
sure time of 4.5 h) in an attempt to measure their redshifts. The 300 
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lines mm~^ grism blazed to 5000 A (1.43 A pixel" ^ dispersion) 
and the 400 lines mm~^ grating (1.86 A pixel" ^ dispersion) were 
used with LRIS-B and R, respectively. The 5600 A dichroic was 
used to split the light. A summary of the multiple slit observations 
is given in Table [T] 

The LRIS slit ni asks were reduced using a custom pipeline 
fseelA uger et al.l2008l for details) that bias subtracts, flat-fields, and 
wavelength-calibrates the spectra before extracting and combining 
the one-dimensional spectra for each slit. Redshifts were deter- 
mined for the extracted spectra by employing an automated cross- 
correlation with galaxy template spectra and manually confirming 
or correcting each of these redshifts by eye. In total, 154 indepen- 
dent objects were observed with the slit masks (note that some of 
these objects were also observed previously with long slits). Red- 
shifts were obtained for 84 extragalactic objects and 21 stars were 
identified. 



5 REDSHIFTS 

The combined long and multiple- slit datasets yielded redshifts for 
90 galaxies. The positions, redshifts and HST magnitudes (F555W 
and F814W) of each galaxy are given in Tabled We now give 
a description of the spectra for the two lensing galaxies within the 
Einstein radius of the system and the lensed source. We also discuss 
the properties of the large-scale structure associated with Gl. 



5.1 The lensing galaxies 

Our optical spectrum of the main lensing galaxy Gl is shown in 
Fig. U The spectrum has a spectral shape and several absorption 
features which are consistent with an early-type galaxy with an old 
(red) stellar population. The redshift of Gl, established from eight 
absorption features (identified in Fig.g]), is 0.3648 ± 0.0002. There 
is known to be a powerful radio-loud AGN embedded within Gl 
with radi o-lobes which ext end out to ^ 30 kpc from the core at 
1.4 GHz (iMore et aDl2008h . Interestingly, there is no evidence of 
AGN activity in the form of emission lines, or a strong blue contin- 
uum in the spectrum of Gl. Therefore, we assume that the AGN is 
either obscured at optical wavelengths or is currently going through 
a period where it is optically passive. 

In Fig. [5] we present part of the two dimensional spectrum 
from the slit placed over the second lensing galaxy, G2, during the 
multiple-slit spectroscopic observations on 2005 August 1 (4.5 h 
total integration time). It is clear that the 2-d spectrum is dominated 
by the much more luminous and extended galaxy Gl. There is no 
evidence of any absorption lines over the full wavelength range of 
the spectrum that are at a different redshift from Gl, but we do see 
absorption lines at a similar redshift to Gl in the G2 spectrum (see 
Fig.[5]for an example of Ca II H and K absorption in both spectra). 
This suggests that G2 could be at the same redshift as Gl. How- 
ever, we cannot rule out that there is contamination of the weaker 
G2 spectrum by the spectrum of Gl. Further observations, which 
spatially resolve Gl and G2, will need to be carried out to confirm 
that G2 is a companion to Gl. 



^ The central wavelength was 7974 A for Mask 5, and 8042 A for Masks 6 
and 7. 
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Figure 1. ijop) The gravitational lens system B2 108+2 13 and the position 
of Slit 1 . (Bottom) The slit positions for the long-slit spectroscopy carried 
out across the field of B2108+213 (underlayed is the HST ACS-F814W 
image). Those galaxies in dashed circles are part of the group associated 
with the main lensing galaxy (Gl) at a redshift of 0.365. The area of sky 
shown is 180 X 180 kpc^ at the redshift of Gl. 



5.2 The lensed source 

In Fig. [6] we present the quasar spectrum of lensed image A. The 
blue-side of the spectrum is featureless and is slightly rising from 
^3800 to 5800 A. This spectral shape is consistent with a BL Lac 
type radio source (e.g., see Marcha et al. 1996), whose optical out- 
put is dominated by the non-thermal emission from the beamed 
quasar. The flat radio spectra of the lensed images of B2 108-1-2 13, 
and the lack of extended structure observed in the radio maps 
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Table 2. Positions, redshifts and spectral properties of galaxies with confirmed redshifts. 



Object 




RA 
(J2000) 


Dec 
(J2000) 


Redshift 


Mask 


Emission line(s) 


F555W 


F814W 


J21 1042.75+213246.3 


21 


10 42.750 


21 32 46.27 


0.7033 


2108m2 


absorption lines only 






J21 1045.30+213204.3 


21 


10 45.304 


21 32 04.32 


0.2875 


2108m7 


[Oii],H/3,[Oiii] 


20.91 ±0.03 


19.54 ±0.03 


J21 1046.38+213403.5 


21 


10 46.376 


21 34 03.47 


0.3639 


2108m7 


[On] 






J21 1047.79+212856.4 


21 


10 47.787 


21 28 56.44 


0.3930 


2108m7 


[Oii],H/3 






J21 1047.79+213345.8 


21 


10 47.789 


21 33 45.78 


0.3628 


2108m7 


[Oii],H^,[Oiii] 






J21 1047.80+213137.2 


21 


10 47.797 


21 31 37.23 


0.5281 


2108m7 


[On] 


23.61 ±0.04 


21.10±0.03 


J21 1048.97+213151.9 


21 


10 48.973 


21 31 51.89 


0.7028 


2108m2 


[On] 


24.09 ±0.04 


21.31 ±0.03 


J21 1049.73+212803.2 


21 


10 49.729 


21 28 03.17 


0.5295 


2108m7 


absorption lines only 






J21 1049.80+213341.8 


21 


10 49.795 


21 33 41.76 


0.1405 


2108m7 


[Oii],H^,[Oiii],Ha,[Nii],[Sii] 






J21 1049.84+213053.9 


21 


10 49.839 


21 30 53.85 


0.2727 


2108m7 


[Oii],H7,H/3,[Oiii],Ha,[Nii] 


21.97 ±0.03 


20.38 ±0.03 


J21 1049.99+213318.5 


21 


10 49.992 


21 33 18.53 


0.6539 


2108m7 


[Oii],H/3,[Oiii] 


23.94 ±0.04 


21.81 ±0.03 


J211050.31+213110.1 


21 


10 50.314 


21 31 10.11 


0.3631 


2108m7 


absorption lines only 


24.06 ± 0.04 


22.14 ±0.03 


J21 1050.34+213025.6 


21 


10 50.341 


21 30 25.61 


0.5607 


2108m5 


[Oil] 


23.80 ±0.04 


21.31 ±0.03 


J21 1050.69+213420.6 


21 


10 50.693 


21 34 20.65 


0.3630 


2108m5 


[Oii],H/3,[Oiii] 






J21 1050.86+212948.6 


21 


10 50.856 


21 29 48.62 


0.3667 


2108m5 


absorption lines only 


23.80 ±0.04 


21.11 ±0.03 


J211051.32+212949.1 


21 


10 51.321 


21 29 49.12 


0.3646 


2108m7 


absorption lines only 


23.97 ±0.04 


21.99 ±0.03 


J211051.58+213003.6 


21 


10 51.582 


21 30 03.58 


0.3668 


LSlit 


absorption lines only 


23.37 ±0.04 


21.13±0.03 


J211051.77+212929.1 


21 


10 51.768 


21 29 29.14 


0.6971 


2108m7 


[Oii],H^ 






J211051.83+212936.9 


21 


10 51.826 


21 29 36.94 


0.3656 


2108m5 


[On] 


22.42 ± 0.03 


20.29 ±0.03 


J21 1052.02+213358.2 


21 


10 52.017 


21 33 58.24 


0.3652 


2108m2 


absorption lines only 






J21 1052.44+213236.6 


21 


10 52.438 


21 32 36.57 


0.3674 


2108m6 


absorption lines only 


22.49 ±0.03 


20.16 ±0.03 


J21 1052.52+212722.4 


21 


10 52.515 


21 27 22.36 


0.8197 


2108m7 


[Oil] 






J21 1052.86+212851.2 


21 


10 52.861 


21 28 51.25 


0.3641 


2108m5 


absorption lines only 






J21 1053.25+213055.0 


21 


10 53.250 


21 30 55.05 


0.3657 


LSlit 


absorption lines only 


23.12 ±0.03 


20.84 ±0.03 


J21 1053.41+212931.5 


21 


10 53.414 


21 29 31.45 


0.3640 


2108m6 


[Oii],H/3,[Oiii] 






J21 1053.65+213049.8 


21 


10 53.650 


21 30 49.80 


0.3678 


LSlit 


[Oil] 


22. 14 ±0.03 


19.95 ±0.03 


J21 1053.71+213107.0 


21 


10 53.709 


21 31 07.02 


0.3088 


LSlit 


[Oii],H/3,[Oiii] 


23.68 ±0.04 


22.38 ±0.03 


J21 1053.76+2131 10.2 


21 


10 53.760 


21 31 10.21 


0.3603 


LSlit 


absorption lines only 


21.63 ±0.03 


19.38 ±0.03 


J21 1053.79+213029.5 


21 


10 53.790 


21 30 29.52 


0.3580 


2108m5 


absorption lines only 


22.52 ±0.03 


20.23 ±0.03 


J21 1053.80+213024.7 


21 


10 53.799 


21 30 24.72 


0.3647 


2108m7 


[On] 


22.12 ±0.03 


20.08 ±0.03 


J21 1053.92+212907.1 


21 


10 53.917 


21 29 07.06 


0.3658 


2108m6 


[Oii],H^,[Oiii] 






J21 1054.02+212847.6 


21 


10 54.025 


21 28 47.61 


0.2436 


2108m2 


[Oii],H/3,[Oiii] 






J21 1054.1 1+213057.7 


21 


1054.111 


21 30 57.73 


0.3643 


LSlit 


absorption lines only 






J21 1054.16+212810.4 


21 


10 54.163 


21 28 10.37 


0.1837 


2108m5 


[Oii],H/3,Ha,[Nii],[Sii] 






J21 1054.36+213016.5 


21 


10 54.356 


21 30 16.51 


0.3649 


2108m5 


absorption lines only 


22.35 ±0.03 


20.02 ± 0.03 


J21 1054.42+212936.1 


21 


10 54.418 


21 29 36.13 


0.3917 


2108m5 


[Oii],H/3,[Oiii] 






J21 1054.61+213228.9 


21 


10 54.610 


21 32 28.86 


0.7982 


2108m7 


[Oii],H(5,H7,H/3,[Oiii] 


23.83 ±0.03 


22.23 ± 0.03 


J21 1054.81+212929.5 


21 


10 54.809 


21 29 29.52 


0.3670 


2108m5 


[Oii],H/3 






J21 1054.97+212804.2 


21 


10 54.969 


21 28 04.25 


0.3907 


2108m5 


absorption lines only 






J21 1054.98+213032.5 


21 


10 54.980 


21 30 32.55 


0.3652 


LSlit 


absorption lines only 


22.85 ±0.04 


20.75 ±0.03 


J21 1054.98+212947.2 


21 


10 54.984 


21 29 47.20 


0.4261 


2108m2 


[On] 






J21 1055.03+213102.4 


21 


10 55.034 


21 31 02.41 


0.3581 


LSlit 


absorption lines only 


22.80 ±0.03 


20.62 ±0.03 


J211055. 15+213115.1 


21 


10 55.152 


21 31 15.12 


0.3611 


2108m2 


[On] 


23.52 ±0.03 


21.85 ±0.03 


J21 1055.34+212804.5 


21 


10 55.340 


21 28 04.53 


0.3910 


2108ml 


[Oii],H/3,[Oiii] 






J21 1055.71+213127.7 


21 


10 55.708 


21 31 27.71 


0.3618 


2108m5 


[Oii],H/3 


22.52 ±0.03 


20.35 ±0.03 


J21 1055.82+213352.5 


21 


10 55.815 


21 33 52.52 


0.3600 


2108m2 


[On] 






J21 1055.84+212926.3 


21 


10 55.843 


21 29 26.32 


0.6099 


2108m4 


absorption lines only 






J21 1055.84+213327.0 


21 


10 55.844 


21 33 27.00 


0.3612 


2108m6 


[Oii],[Oiii] 






J21 1055.92+212921.7 


21 


10 55.922 


21 29 21.70 


0.8485 


2108m5 


[On] 






J21 1056.00+213347.3 


21 


10 55.997 


21 33 47.34 


0.3609 


2108m5 


[On] 






J21 1056.02+212803.1 


21 


10 56.019 


21 28 03.11 


0.3533 


2108m7 


[Oii],H^,[Oiii] 






J21 1056.09+212816.4 


21 


10 56.092 


21 28 16.42 


0.3652 


2108m7 


[Oii],H/3,[Oiii] 






J21 1056.13+213432.8 


21 


10 56.133 


21 34 32.77 


0.2991 


2108m5 


absorption lines only 






J21 1056.35+213413.5 


21 


10 56.354 


21 34 13.53 


0.7962 


2108m6 


[On] 






J21 1056.40+213218.6 


21 


10 56.400 


21 32 18.65 


0.3630 


2108m6 


[On] 


20.74 ±0.03 


18.73 ±0.03 


J21 1056.62+213401.9 


21 


10 56.618 


21 34 01.93 


0.3625 


2108m4 


absorption lines only 






J21 1056.82+213158.9 


21 


10 56.824 


21 31 58.92 


0.3634 


2108m5 


absorption lines only 


22.88 ±0.03 


20.60 ±0.03 


J21 1056.86+213452.7 


21 


10 56.862 


21 34 52.73 


0.3639 


2108m6 


[On] 






J21 1056.98+212918.5 


21 


10 56.978 


21 29 18.54 


0.3925 


2108m7 


[Oii],H/3,[Oiii] 






J21 1057.14+212749.4 


21 


10 57.137 


21 27 49.41 


0.6838 


2108m6 


[Oii],H/5 






J21 1057.47+212759.8 


21 


10 57.473 


21 27 59.80 


0.3537 


2108ml 


[Oii],[Oiii] 






J21 1057.50+213349.4 


21 


10 57.501 


21 33 49.44 


0.3628 


2108m6 


[On] 






J21 1057.52+213335.0 


21 


10 57.521 


21 33 35.04 


0.3619 


2108m5 


[On] 
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Table 2 - continued 



Object 




RA 
(J2000) 


Dec 
(J2000) 


Redshift 


Mask 


Emission line(s) 


F555W 


F814W 




21 


10 57.547 


91 1 54 

Zi JJ iU.JH- 


0.3651 


91 08m^S 

ZiL/OiiiU 


ro iii WR Wrv 






T9 1 1 057 58-1-9 1 '^409 7 


21 


10 57.575 


91 '^4 09 71 

Z i JH- L/Z. / i 


0.3657 


91 08m^S 

ZiL/OiiiU 


ro Til 
\SJ lij 






T9 1 1 057 1 -1-9 1 ^ 1 


21 


1 57 ^S1 1 

iU J / .Ui i 


91 ^1 5^S 04 
z i J i J u . yj'-T 


0.3680 


91 08m^S 

ZiL/OiiiU 


aUSUipiiUil iiilCS Uiliy 




91 1 n -1- 0'^ 

Zi . lU ZIZ U.UJ 


T91 1 057 ^Szl-i-91 ^955 8 


21 


1 57 fi^'> 


91 ^9 55 lf\ 
Z i JZ J J . / u 


0.3566 


91 08m 5 

ZiL/OiiiJ 


aUSUipiiUil iiilCS Uiliy 






T9 1 1 057 80-1-9 1 ^ 1 ^Zl 


21 


1 57 80'^ 
iu J / .oKjj 


91 ^1 ^4 ^^1 

Zi Ji J^-.Oi 


51 95 

U.J iZJ 


9 1 08m^^ 
ziuoiiio 


ro ill wF\ rOiiii 

\\j iij,np,Lw iiij 


90 no _|_ n 04 
AD.yo ZC U.U'4- 


99 m _|_ n 0'^ 

ZZ.UJ ZIZ U.UJ 


T9 1 1 057 89-1-9 1 ^OzlQ 9 


21 


1 57 89'^ 

iU J / .OZJ) 


91 '^0 4Q 99 

Z i JU H-7.ZZ 


0.3629 


91 08m7 

ZiL/Oiii / 


rOiii wR rOiiii nviii r<^iiii 

\\j iij ,n/J , \\J iii J ,nci, iij , iij 


9a r)9 -1- 0^ 

ZJ.UZ ZC U.UJ 


91 _|_ n 0^ 

Zl.H-U ZIZ U.UJ 


T9 1 1 057 8^-1-9 1 ^Ozl^ 7 


21 


1 57 897 
iu J / .oz / 


91 ^0 4^ ^^Q 
Zi jyj ^j.yjy 


8959 
u.ozjz 


91 08m^^ 
ziuoiiio 


Hbsorption lines only 


94 ^^8 -1- 0^^ 


91 -1- 0^ 

Zi . JU ZIZ U.UJ 


T9 1 1 058 00-1-9 1 ^ 1 44 Q 


21 


10 58 00^ 


91 '^1 44 8Q 
Zi ji '-t'-t, oy 


'^588 

L/. J JOO 


9 108m 5 

Z i VJOliiJ 


clUaUipLiUll iillCS Uiliy 


99 Q1 _L A AO 
ZZ.!7i ZIZ U.UJ 


90 _L A AO 
ZU.UU ZIZ U.UJ 


T9 1 1 058 08-1-9 1 ^^05 Q 


21 


1 58 081 
iu Jo.uoi 


91 ^^05 Q9 
Zi DD yjJ.yz. 


'^585 
u. J JoJ 


91 08m 5 
ziuoiiij 


dUsUipLiuii lilies uiiiy 






T91 1058 5^^-1-9 19Q09 8 

JZi iUJO. JU+ZiZ^/L/Z.O 


21 


1 58 558 


91 9Q 09 lf\ 

Z i Z7 L/Z. / U 


0.3643 


91 08m^S 

ZiL/OiiiU 


dusuipLiuii lilies uiiiy 






T9 1 1 058 5^^-1-9 1 ^0^5 

JZi iUJO. JO+Zi J>UJ>J.U 


21 


1 58 5^^0 
iu Jo.jou 


91 ^0 ^5 01 

Zi DkJ JJ.Ui 


U. JOOZ 


91 08m^^ 
ziuoiiio 


absorption lines only 






T91 1 058 71 -1-91 ^^9^S Q 

JZi iUJO. / i+Zi JJZU." 


21 


1 58 71 5 
iU Jo. / i J 


91 9^S Q4 

Zi JJ ZU.7H- 


0.1682 


91 08m4 

ZiL/OiiiH- 


\\j iij,nci,Li> iij,LO iij 






T91 1 058 Q9-1-91 ^941 

JZi iUJO."Z+Zi J)ZH-i .U 


21 


1 58 Q1 8 

iU Jo." i o 


91 ^9 40 Q7 
Z i JZ '-tKJ.y 1 


0.3607 


91 08m5 

ZiL/OiiiJ 


dUSUipLlUll lilies Uiliy 






T91 1 05Q 1 0-1-91 ^951 1 

JZi LKJJy. iL/T^Zi J)ZJ i . i 


21 


1 5Q OQ^S 


91 ^9 51 15 

Zi JZ J i . i J 


0.3591 


91 08m5 

ZiL/OiiiJ 


dusuipLiuii lilies Uiliy 






T9 1 1 05Q ^^0-1-9 1 ^01 Q 1 


21 


1 5Q ^^01 

iU J^.OUi 


91 ^0 1 Q 9Q 

Zi DKj Vy.ZJy 


1 481 5 

i .'4-0 i J 


91 08m^^ 
ZiL/OliiO 


absorption lines only 






T9 1 1 05Q 89-1-9 1 ^ 1 04 7 

JZi iUJ".OZT^Zi J) iUH-. / 


21 


1 5Q 81 

iU J!7.0 iU 


91 '^1 04 74 

Z i J i L/H- . / H- 


0.3649 


91 08m^S 

ZiL/OiiiU 


dusuipLiuii lilies Uiliy 


99 04 _|_ n m 

ZZ.OH- ZIZ U.UJ 


9A 4A _|_ A Aa 
ZU.H-U ZIZ U.UJ 


T9 1 1 05Q 8^^-1-9 1 ^9^^^ 1 

JZi iUJ^.OO+Zi J>ZJ>0. i 


21 


1 5Q 8^^5 

iU J7.0OJ 


91 '^9 08 
Zi JZ jo.uo 


^5QQ 


91 08m 5 
ziuoiiij 


absorption lines only 






T91 1 100 08-1-9 n954 1 

J Z i i i VJVJ.VJOT^Z i JZ Jt". i 


21 


1 1 00 085 

i i VJL/.L/OJ 


91 ^9 54 1 ^ 

Zi JZ jT-.iJ 


0.3168 


91 08m^S 

Z i VJOliiU 


ro Til 
\\j iij 






T9 1 1 1 00 ^S7-i-9 1 ^ 1 57 Q 

JZi i iVJL/.U / T^Zi J i J / .7 


21 


1 1 00 674 


91 '^1 57 Q'^ 
z i J i J / .y J 


0.3639 


91 08mfS 

Z i VJOliiU 


ro Til 
\\j iij 






J211100.73+212919.8 


21 


11 00.732 


21 29 19.76 


0.5123 


2108m6 


[On] 






J21 1100.86+212838.0 


21 


11 00.859 


21 28 38.04 


0.6456 


2108m6 


[Oii],H^,[Oiii] 






J21 1100.98+212825.6 


21 


11 00.975 


21 28 25.64 


0.8275 


2108m6 


[Oil] 






J21 1100.98+213249.7 


21 


11 00.984 


21 32 49.68 


0.0837 


2108m4 


Ha,[S II] 






J21 1103.21+213020.7 


21 


11 03.213 


21 30 20.68 


0.5880 


2108ml 


[Oil] 






J21 1103.33+212938.3 


21 


11 03.331 


21 29 38.34 


0.3665 


2108m4 


[Oii],H^,[Oiii] 






J211104.15+212856.7 


21 


11 04.150 


21 28 56.72 


0.2314 


2108m2 


[O II],H7,[0 iii],Ha,[Nii],[S II] 






J211104.21+212912.4 


21 


11 04.206 


21 29 12.42 


0.2325 


2108m6 


[Oii],H/3,[Oiii],Ha,[Nii],[Sii] 








4000 4500 5000 5500 6000 6500 7000 7500 8000 

Wavelength (A) 

Figure 4. The optical spectrum of the main lensing galaxy (Gl) of B2 108+2 13 taken with LRIS on the W. M. Keck Telescope on 2003 October 24. Gl is an 
early-type galaxy at a redshift of 0.3648 ± 0.0002. The main absorption features used to determine the redshift of Gl have been identified and the skylines 
have been labeled with the crossed circles. 



dlVlcKean et alJl2005l :MQreetal are both in agreement with 

the classification as a beamed radio source viewed down the axis 
of the radio-jet. The red side of the spectrum shows a hint of emis- 
sion from the quasar host galaxy, but also has contamination from 
the wings of lensing galaxy Gl (see Fig.©. In Fig.|7]we show the 
spectrum of the extended arc of the host galaxy of the lensed quasar 
- note that this spectrum has not been flux-calibrated, so the data 
are shown in relative counts without the shape of the bandpass be- 



ing corrected for. There is evidence of a break in the spectrum at 
-6700 A, which we identify as the B aimer break. Note that the 
break in the spectrum cannot be the Lyman break since there is no 
evidence of absorption in the quasar spectrum below 6700 A by 
the Lyman forest. We also see a step in the spectrum of the lensed 
quasar at 6660 A in the response corrected spectrum of image A 
shown in Fig. [6] Consistent with the B aimer break in the spectrum 
is the tentative detection of the Mg-b absorption line at 8661 =t 27 A 
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Figure 2. The V-/? colour band magnitude plot of the B2 108+2 13 field. 
The red diamonds mark the observed colours of Gl and the 4 closest galax- 
ies to Gl detected from the long-slit observations, whereas the blue circles 
mark those objects with measured redshifts. Note that this relative colour 
plot is made from non-photometric data and the expected V—R colour of an 
elliptical galaxy at redshift 0.365 is 1.2. 
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Figure 3. The R-I colour / band magnitude plot of the B2108+213 field. 
The red diamonds mark the observed colours of Gl and the 4 closest galax- 
ies to Gl detected from the long-slit data, whereas the blue circles mark 
those objects with measured redshifts. Note that this relative colour plot is 
made from non-photometric data and the expected R—I colour of an ellip- 
tical galaxy at redshift 0.365 is 0.9. 
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Figure 5. The two dimensional spectrum of the B2 108+2 13 lens system 
between 5330-5450 A. The position of lensed image A is marked as the 
origin of the spatial axis. The spectrum shows the clear Ca ii H and K ab- 
sorption lines in lensing galaxy Gl (— 3.2 arcsec from lensed image A). The 
spectrum at the position of the second lensing galaxy (—2.2 arcsec from im- 
age A) shows a hint of the same Can H and K absorption lines. However, 
the signal-to-noise ratio and spatial resolution of the spectrum is not good 
enough to confirm that G2 is at the same redshift as Gl. 
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Figure 7. The optical spectrum (5700-9200 A) of the gravitational arc from 
the lensed image A. The spectrum shows a possible spectral break and Mg-b 
absorption, giving a likely source redshift of 0.674 ±0.005. The spectrum 
has been smoothed with a box-car of width 20 A and has not been flux- 
calibrated or corrected for the spectral response of the instrument. 



in Fig.|7] This gives a likely redshift of Zs = 0.674 ± 0.005 for the 
lensed source. Interestingly, there is an absorption line seen in the 
lensed quasar spectrum at 7172 zb 17 A which could correspond 
to the G-band at redshift 0.666 ± 0.004. However, this line is very 
close to a deep absorption line seen in the lens galaxy spectrum, 
which could be contaminating the absorption seen in the lensed 
quasar spectrum. There is no evidence of absorption line systems 
in any part of the spectrum of the lensed quasar, which would agree 
with the source redshift being relatively low. 



5.3 Large-scale structure velocity dispersion 

The redshift distribution of the 90 extragalactic objects found from 
the spectroscopic observations is shown in Fig. [8] There is a clear 
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Figure 6. The optical spectrum of the lensed quasar (image A) of B2 108+2 13 taken with LRIS on the W. M. Keck Telescope on 2003 October 24. The flat 
optical spectrum of image A is consistent with a BL Lac type radio source. The red end of the spectrum shows some emission from the quasar host galaxy and 
contamination from the lens galaxy, Gl. Sky emission and absorption lines have been labeled with the crossed circles. The spectrum has been smoothed using 
a boxcar of 7 A. 



spike at redshift 0.365, the redshift of the lensing galaxy Gl. The 
galaxy membership of the large-scale structure was found by ap- 
plying a linking velocity kernel. We define the linking velocity as 
an initial guess at the velocity dispersion of the structure. This 
linking velocity is also used to estimate the structure member- 
ship; the candidate galaxies are ordered by velocity and each mem- 
ber must be within the linking velocity of the next member. The 
group n iembership is refine d usi ng the iterat i ve pr ocedure out- 
lined in IWilman etaP (l2005h and lAuger etaP ilOOH) . Due to the 
moderately large number of structure members, the bi- weight esti- 
mate of scale is used to determine the effecti ve group velocity dis- 
persion (see Beer s. Flvnn & GebhardtlfToOOl) . Several linking ve- 
locities were used to estimate the group membership (450, 650 
and 1100 kms~^; see Table [5]), which reaches a maximum of 52 
members using a linking velocity of 1100 kms~^. The three link- 
ing velocities returned group velocity dispersions of 451, 489 and 
694 km s~^, respectively, and in each case the velocity distribution 
is non-Gaussian. This suggests that the structure is not viralised and 
is still in the process of forming. 

In Fig. 121 we show the positions of the 52 galaxies identified 
using the 1100 kms~^ linking kernel. Also, the colour of each 
galaxy indicates how much the galaxy is blue or redshifted with 
respect to the systemic velocity of the structure. Interestingly, there 
appears to be an over-density of blueshifted galaxies to the north- 
east of the lens galaxy, and an over-density of redshifted galax- 
ies to the southwest. Also, the divide between these blue and red- 
shifted compo nents is straddled by the X-ray emission detected 
with Chandra Jpassnacht et al.ll2008l) . There is also no evidence of 
X-ray emission across the spatial distribution of the spectroscop- 
ically confirmed members of the lensing structure. It may be that 
we are observing the ongoing merger of two independent struc- 
tures in this system. Based on this conclusion, we have identified 
two separate groups of galaxies; one which is north of the X-ray 
emission, and the other to the south. The velocity distributions of 
these two groups are shown in Fig.[TO]as the blue and red distribu- 
tions. Further observations will certainly have to be carried out to 
investigate the possibility that there are two merging groups in this 
system. For example, our redshift data at present are biased along 
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Figure 8. The galaxy redshift distribution from the long and multi-slit 
datasets of the 82108+213 field. Redshifts were obtained for 90 extragalac- 
tic objects. 



an almost north- south orientation, with the lens system at the cen- 
tre. This is because our imaging data with LRIS are orientated in 
this direction (see Table[T]). 



6 STELLAR KINEMATICS AND MORPHOLOGY OF Gl 

The long- slit spectrum of Gl was taken with a 1 arc sec wide slit 
over an integration time of 1.5 h. This set-up produced a spectrum 
with a high enough spectral resolution (^174 km s~^ FWHM) and 
signal-to-noise ratio (^^80) for the line-of- sight stellar velocity dis- 
persion of Gl to be measured. The stellar velocity dispersion pro- 
vides an additional constraint to the mass model and allows the 
average logarithmic matter density profile of the lensing galaxy to 
be measured when combined with the mass estimate from grav- 
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Figure 9. The position of the group galaxies relative to the group centre for 
the 1100 kms~^ linking kernel. The position of the lensing galaxy (Gl) 
is shown by the cross. The intensity of the blue/red colour of each galaxy 
is dependent on the velocity relative to the systemic velocity of the group. 
The conto urs show the distribution of the X-ray emission around the lensing 
group [see lFassnacht et al] i2008h for details] and the circle marks a radius 
of715kpc. 
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Figure 10. The velocity distribution for the B2 108+2 13 large-scale struc- 
ture using a linking velocity kernel of 1100 kms~^ (white). The velocity 
distributions of the blue and redshifted structures is also shown. 



itational lensing (e.g IXreu & Koopmansl l2004l : ISand et all [2004 
iKoopmans et alJl2006h . 

The stellar velocity dispersion of Gl w as measured using tem- 
plate spectra of eight G and K type stars dSand et aljbooi) . The 
templates were first smoothed to the resolution of the Gl spectrum, 
before the best-fitting velocity dispersion was determined from a 
linear combination of a polynomial continuum and the ensemble 
of templates. The measured stellar velocity dispersion within an 
aperture of 0.5 arcsec diameter is (Jap = 325 =b 25 kms~^. This is 



Table 3. The properties of the galaxy group associated with B2 108+2 13 
using a kernel of width 450, 650 and 1 100 km s~ ^ . The right ascension and 
declination of the galaxy group are calculated from the mean position of the 
group galaxies. 



Kernel 


RA 


Dec 


Number of 






(J2000) 


(J2000) 


galaxies 


(kms-1) 


450 


21 10 55.147 


21 31 36.92 


41 


451 ±45 


650 


21 10 55.108 


21 31 37.64 


44 


489 ±47 


1100 


21 10 55.382 


21 31 32.09 


52 


694 ± 93 



much larger than the characteristic velocity dispersion for the early- 
type lens galaxy population detern iined from gravit ational lensing 
statistics (a* = IQStg^ km s' MChae et aDl2002h and fr om the 
velocity dispersion function of galaxies ('Sheth et aL""2003"). How- 
ever, the large value is typical of the most massive early-type galax- 
ies that are usually f ound at the centres of groups and clusters (e.g. 
iBernardi et aDl2007h . 

The analysis of the high resolution optical and infrared imag- 
ing of B2 108+2 13 with the HST found that the surface brightness 
dist ribution of Gl was w ell represented by a de Vaucouleurs pro- 
file (iMcKean et al.ll2005h . However, fitting a single de Vaucouleurs 
profile to galaxy Gl revealed asymmetric residuals that suggested 
a more complex model for the light distribution was required. The 
data presented here have shown Gl to be part of a group or cluster, 
that it is the brightest group/cluster galaxy and has a large cen- 
tral stellar velocity dispersion. Also, Chandra imaging showed that 
Gl is close to the centre of the exten ded X-ray emission from the 
group/cluster jpassnacht et al.ll2008[). Thus, it is like lv that Gl is 
a giant elliptical cD galaxy fe.g. lMatthews et al ][l964). Therefore, 
we have reanalysed the surface brightness profile of Gl, using a 
new model that takes into account a large central bulge and an ex- 
tended stellar envelope, typical of cD galaxies. This reanalysis is 
required because the effective radius of the bulge surface brightness 
profile is used in the dynamical models for the mass distribution in 
Section |7] 

The surface brightness profile fitting of the lensing galaxy 
light from the ACS F814 W and NICMOS F160W imaging was car- 
ried out using GALFIT ( |Peng et al. 2002). The bulge of lensing 
galaxy Gl, and all of the light from lensing galaxy G2 were fitted 
using Sersic profiles, with the index set to n = 4 (de Vaucouleur 
models). The extended stellar envelope of Gl was fitted using a 
Sersic profile, with the index left as a free parameter. The profiles 
were convol ved with the HST point spread function prior to fitting 
(lKristlll995h . The extended arc emission from the lensed images 
was masked out. 

Including a stellar envelope component results in a signifi- 
cantly better fit to the observed surface brightness distribution. For 
the ACS F814W dataset, the reduced of the model with a single 
de Vaucouleurs profile fitted to Gl is 1.8, whereas also including 
an envelope in the model lowers the to 1.3. Note that fitting a 
two component model to Gl also removed the asymmetric resid- 
uals that were seen a round the galaxy in the previous analysis by 
iMcKean et al ] (l2005h . The fitted Sersic slope of the stellar enve- 
lope is 1.4 and the position is offset from the stellar bulge of Gl by 
^1.1 arcsec to the east. These results are both consistent with the 
morphology of envelopes in cD galaxies, which are often flattened 
and offset from the bulge. The parameters of the fitted light distri- 
butions to the F814W data are given in Table |4] Similar results are 
found for the NICMOS F160W dataset. However, the values of the 
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Table 4. The surface brightness distributions of the two galaxies within the 
Einstein radius (Gl and G2). Only the parameters for the F814W imaging 
are presented because the field of view of the ACS is large enough to include 
all of the fight from the extended stefiar envelope of Gl. 



Parameter 




Gl 


G2 


Unit 




bulge 


envelope 






Total magnitude 


17.74 


17.95 


20.42 


Vega-mag 


Effective radius 


1.17 


6.56 


0.24 


arcsec 


Slope 


4.00 


1.44 


4.00 




Axial ratio 


0.88 


0.81 


0.77 




Position angle 


61 


161 


128 


degrees 



slope (0.7) and effective radius (2.1 arcsec) for the stellar envelope 
were significantly different because the field of view of NICMOS is 
too small to properly fit the extended envelope component of Gl. 
We find that by using a two component model for Gl, the effective 
radius of the bulge is decreased by 16 per cent, whe n compared to 
what was found from the single component model (iMcKean et alj 
l2005h . 



7 GRAVITATIONAL LENS MASS MODEL 

In lMore et alJ (l2008h it was shown that a simple mass model con- 
sisting of two singular isothermal spheres (SIS) plus a moderate 
external shear (7 ^ 0.04) could fit the observed VLBI constraints 
from the B2 108+2 13 lensed images reasonably well, without the 
need for a large group or cluster contribution to the local shear 
and/or convergence. Despite this success, the flux-ratio of the two 
lensed images A and B could only be well-fit if the density pro- 
file of Gl (and G2 to a lesser extent) was steepened to a loga- 
rithmic density slope of 7^ = 2.45lQ }g (68 per cent confidence 
level). This is relatively steep, but not inc onsistent with some of 
the sy stems found in the SLAGS survey jKoopmans et aljEood 
12009) or, for example, with PG 1115-1-080, a nother l ens system 
close l y associated with a co mpact group (e.g. iTreu & Koopman s 
l2002l : lMomcheva et al.ll2006h . 

Here, we first attempt to include the contribution of the bright- 
est nearby group members of B2 108-1-2 13 into the gravitational 
lensing mass model in order to assess their effect, before incor- 
porating the stellar dynamics of Gl and the extended arc structure. 
We also assess the effect of including a mass-sheet to account for 
the group. In addition, we attempt more general models with shear 
and a constant convergence gradient {dn/dr) over the extent of the 
lensed images. The combination of these two can closely mimic 
the general effects of an external group or cluster for which a more 
precise mass model has yet to be established. To assess the relative 
effects of Gl and G2, and that of the group/cluster, we study two 
classes of mass models, namely, (i) isothermal (7^ = 2) and (ii) 
power-law density profiles; in both cases we investigate spherical 
and elliptical models. 

7.1 Models including the population of sub-haloes 

We include the field galaxies as follows. For both the isothermal 
and non-isothermal cases, we normalize the masses inside the Ein- 
stein radius of each group member (including G2) to that of Gl, 
using the proportionality relation M oc L^^^, where M is the mass 
within a fixed aperture for an isothermal mass model. This case 
is consistent with the Faber-Jackson relation jFaber & JacksonI 



[1973). 

All galaxies (besides G2) are assumed to be spherical to 
a first approximation. We also only include in our models those 
group members within 45 arcsec of the lens system, since out to 
this distance, we have good photometry and positions for the group 
galaxies from the HST ACS data. This adds eleven galaxies to the 
mass models. Even though there are galaxies found farther from the 
lens system, their contribution to the overall convergence is consid- 
ered to be neglible; recall that the convergence n ^ 6£;/r for an 
isothermal mass model, where r is the distance from the galaxy 
and is the Einst ein radius. We use exactly the same constraints 
on the models as in iMore etaLl(l2008h . As a sanity check, we have 
tested o ur code against their results [for which the code of iKeetonI 
(l2001ah was used] and confirm all their findings and values to 
within the expected numerical errors. 

Adding the field contribution to the models in the form of 
independent sub-haloes, besides an overall lowering of the lens- 
strength of both Gl and G2 (see Table[5]) and an increase of the ex- 
ternal shear, does not alter any of the model parameters outside the 
error range. The isothermal model with shear (2SIS-i-shear-i-field) 
h as a reduced = 7.0, which is almost identical to the model 
of lMore etaD ( l2008h where the field contribution was not included 
(reduced — 6.9). Like before, the model fits the positions of the 
lensed images A and B well, but fails to recover the observed flux- 
ratio. Therefore, solely adding the contribution of the field galax- 
ies as sub-haloes fails to reproduce the observed properties of the 
lensed images. The additional convergence of the field results in a 
lowering of the Einstein radii of Gl and G2 by about 15 per cent. 

We find that the positions and flux-ratios of the lensed 
images can both be reproduced if we invoke either a con- 
stant convergence gradient of dn/dr — 0.059 ±0.002 arcsec"^ 
(2SIS-i-shear-i-field-i-gradient model), or a steeper density profile for 
the lensing galaxy Gl of 7' = 2.51 ±0.01 (2PLS-hshear-hfield 
model) - resulting in a reduced (see Table [5]). However, 

it was not possible to differentiate between these two models based 
on the constraints provided by the VLBI observations of the lensed 
images. In Fig. [TT] we show the critical curves and caustics for 
the model which includes an external shear, the field and a con- 
stant convergence gradient - the critical curves for the other two 
models investigated here are similar. Note that the constant con- 
vergence gradient has a position angle of 15 degrees east of north 
(i.e. north-south), which is approximately increasing in the direc- 
tion toward the centre of the X-ray emission found from this system 
(iFassnacht et alJl2008h . 



7.2 Models using lensing and stellar dynamics 

We now include in our analysis the stellar kinematic data to de- 
termine the average logarithmic density slope between the Ein- 
stein and effective radii of Gl. The inputs to the kinematic model 
are the estimated stellar velocity dispersion (325 ±25 kms~^) 
and the effective radius (1.17 arcsec) of Gl. These models are 
constrained using the mass within the Einstein radius obtained 
from the gravitational lensing analysis. For the Einstein radius 
of Gl, we have adopted a value of bE = 1.51 ±0.26 arcsec, 
which was taken from the 2SIS-i-shear-i-field model (see Table [5]). 



^ The convergence gradient is defined as di^/dr and is added to the model 
at the position of the lensing galaxy, where its absolute value is 0. The 
strength of the gradient and its direction add two free parameters to the 
model. 
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Figure 11. The critical curves (red) and caustics (blue) of the mass model 
for B2 108+2 13 that includes the field galaxies, an external shear and a con- 
vergence gradient. 



Note that we have increased the uncertainty on the Einstein ra- 
dius of Gl to account for the possibility that G2 is more or less 
massive than has been assumed using the Faber- Jackson relation. 
The projected masfl within the Einstein radius of Gl is found 
to be 7.2 ±2.5 x 10^^ Mq. Using this mass and its uncertainty, 
the spherical Jeans equati ons are solved fo llowing the method de- 
scribed in, for example, iKoopmans et al.l 12006). We have also 
tested for the effect of anisotropy in the stellar orbits (/3 = ± 0.25) 
and of using either a Hernquist or Jaffe model for the stellar com- 
ponent of Gl. 

In Fig. [T2I we show the -logarithmic density slope plots 
for lensing galaxy Gl. For each of the combinations of mass, stel- 
lar model and anisotropy, the slope of the density profile is close 
to isothermal. We see that the largest scatter in the derived slope 
comes from the possible range of masses that Gl can have due to 
the uncertainty of including galaxy G2 in the model. The best fit- 
ting density slope is 7^ = 1.981q 5o for the Hernquist model and 
7' = 2.00t 07 for the Jaffe model (both at the 68 per cent confi- 
dence level). The range of slopes allowed for the possible masses 
that Gl can have is 1.82 ^ 7^ ^ 2.16 (68 per cent confidence level). 
Therefore, it is clear from the lensing and dynamics analysis of the 
mass distribution presented here that the average logarithmic den- 
sity slope of Gl is consistent with an isothermal mass model. We 
note that the combined luminous and dark matter density slope of 
7^ = 2.51 required to reproduce both the positions and flux-ratios 
of the lensed images by the lensing model alone would require the 
stellar velocity dispersion of Gl to be of order ^ 600 km s~^. 



^ The mass within the Einstein radius has been calculated assuming an 
isothermal mass distribution; Me = 1.24 x 10^^ 6| (DiDs/Dig) Mq, 
where Di , Dg and Dis are the angular diameter distances to the lens, source 
and between the lens and source, respectively. 
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Figure 12. The best fitting logarithmic density slope of the mass distribution 
for Gl from a joint lensing and dynamics analysis. {Top) The thick black 
lines are for an isotropic distribution of stellar orbits, whereas, the thin grey 
lines correspond to an anisotropy of —0.25 and +0.25. (Bottom) The thick 
black lines are for the mass of Gl assuming the Faber- Jackson relation 
holds for G2. The thin grey lines are for the cases were G2 is more or less 
massive than predicted from the Faber-Jackson relation. (Left and right) 
The analysis for Hernquist (HQ) and Jaffe (JF) models for the stellar mass 
distribution of Gl. 



7.3 Models using the extended arc emission 

We now test models using the extended arc emission from the back- 
ground quasar host galaxy that is seen in both lensed images. In re- 
cent years, sophisticated grid-based modeling algorithms have been 
developed, allowing the properties of the lensing potential to be 
found and for the background object to be rec onstructed from the 
extended source intensity distribut i on (e.^. [ Warren & Dvell2003l : 
iKoopmansI I2OO5I : ISuvu et all [200^ : IVegetti & Koopmansl bOOQah . 
Here, we attempt to constrain the radial profile of the lensing 
mass distribution by investigating the radial structure of the ex- 
tended arc emission. For example, for an isothermal density pro- 
file, the r adial stretch ing shoul d be identica l for both lensed im- 
ages (e.g. IWucknitz. Biggs & Brownd 1200 4V The modelling was 
ca rried out with the adaptive a nd fully Bayesian grid-based code 
of Vegetti & Koopmansl (|2009a). We use for our analysis the NIC- 
MOS F160W data for B2108-h213 that has the light from the lens- 
ing galaxies removed and any residuals masked out. We have also 
masked the bright quasar emission for two reasons. First, the quasar 
emission only supplies limited constraints to the mass model while 
introducing a strong curvature in the reconstructed source. Second, 
in the modelling, we have used a truncated version of the psf, which 
typically produces negligible effects below the noise for most pix- 
els in the image, but would produce non-negligible artifacts for 
those pixels containing the strong point-source quasar emission. 
Ideally we would like to differentiate between the isothermal mod- 
els predicted from the lensing and dynamics analysis and the mod- 
els with a steep power-law slope of 7^ ^ 2.5 found from the VLBI 
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Table 5. The parameters of the gravitational lens models for B2 108+21 3. All positions are given relative to lensed image Al and all position angles are 
measured east of north. The stated model parameters are the median values obtained by carrying out a MCMC simulation and the resulting uncertainties have 
been found from the 68 per cent confidence levels. Note that these are statistical uncertainties and do not reflect any uncertainties due to the choice of models 
used. 



Model 


2SIS-hshear-hfield 


2SIS-hshear-hfield+gradient 


2PLS+shear+field 




Gl position (Aa, A(5) 


1431 ± 1, -2888 i 


:1 1431 ± 1, -2888 ±1 


1431 ± 1, -2888 ±1 


mas 


Gl Einstein radius 


1511±6 


1425 ± 3 


1893 ±11 


mas 


Gl slope 






2.51 ±0.01 




G2 position (Ao;, A(5) 


1273 ±5, -1786i 


:5 1273±4, -1786±4 


1273 ±5, -1786 ±5 


mas 


G2 Einstein Radius 


437 


415 


549 


mas 


External shear 


0.065 ±0.001 


0.040 ±0.001 


0.088 ±0.005 




External shear position angle 


95 ±2 


92 ±2 


77 ± 1 


degrees 


Gradient 




59 ±2 




mas"-*^ 


Gradient position angle 




15±4 




degrees 


Source 1 position (Ao;, A5) 


956, -3542 


1010, -3646 


741, -1957 


mas 


Source 1 flux-density 


1.1 


0.8 


3.8 


mJy 


Source 2 position (Ao;, A(5) 


955, -3542 


1008, -3645 


739, -1956 


mas 


Source 2 flux-density 


0.4 


0.3 


1.3 


mJy 


Source 3 position (Ao;, A(5) 


953, -3543 


1007, -3645 


736, -1958 


mas 


Source 3 flux-density 


0.2 


0.1 


0.6 


mJy 


Flux-ratio (Bl/Al) 


0.63 


0.45 


0.46 




Reduced 


6.96 


0.01 


0.07 





lensing data alone. Therefore, we have used the Bayesian evidence 
values to estimate the best fitting model. 

In Fig. [13] we show the results of fitting an isothermal sphere 
mass model and a spherical power-law model with a density slope 
of 7^ = 2.50. For the case with an isothermal mass model, the Ein- 
stein radii of lensing galaxies Gl and G2 are 1.87 and 0.39 arcsec, 
respectively. Note that for these models the mass-ratio of Gl and 
G2 has been left as a free parameter. The shear strength is found to 
be 0.03 at a position angle of 88 degr east of north. These model pa- 
rameters are similar to, but not formally consistent with the results 
obtained using only the VLBI constraints from the lensed images 
(iMore et alJl20 08). The source has an elliptical shape with the peak 
surface brightness at the centre of the reconstructed host galaxy, 
which is plausible. For the spherical power-law model, the Einstein 
radii of Gl and G2 are lowered to 1.65 and 0.40 arcsec, respec- 
tively, and the shear is increased to 0.075 at a position angle of 
82 degr east of north. Again, the source reconstruction looks plau- 
sible. 

Both the isothermal and power-law models result in a good 
reproduction of the observed source intensity distribution (see 
Fig. [13]), although there are still some systematics left in the im- 
ages. Formally, the power-law model is a better fit to the data, 
with a Bayesian evidence d ifference of A logio E — ^1 [see 
IVegetti & KoopmansI (l2009al) for a discussion on Bayesian evi- 
dence values and gravitational lens modelling]. Note that we did 
not include a convergence gradient for the isothermal case during 
this analysis because here we only wanted to differentiate between 
the two possible density slopes using the extended infrared emis- 
sion from the lensed source. 



7.4 Models including a mass-sheet 

As a final test, we have included a flat mass-sheet into the lens 
model to account for any external convergence that is not taken 
care of through the inclusion of the satellite galaxies of the group. 
The effect of a mass-sheet would be to decrease the amount of mass 
attributed to Gl in the lens models, which would in turn result in 



the logarithmic density slope of Gl to increase for a fixed stellar 
velocity dispersion. This degeneracy between the slope of the den- 
sity profile and the presence of an underlying invisible mass-sheet 
is often referred to as the mass-sheet degeneracy, and can have 
important implications for determining, f or example, th e Hubble 
constant from gravitational lenses (see e.g. lSuvu et a l. 2010). How- 
ever, with the measurement of the stellar velocity dispersion of the 
galaxy, this degeneracy can be broken and the average slope of the 
inner density profile can in principle be determined (see Section 
17.21 ). Alternatively, we can use the slope required from the grav- 
itational lensing model and the measured stellar velocity disper- 
sion to calculate the external convergence from the mass-sheet that 
is needed to reconcile the two. For this model, we use the den- 
sity sl ope 7 ' = 2.45t°:l^ determined from the VLBI constraints 
( More et al ][200 8). with a Gaussian prior on the slope. We also use 
the stellar velocity dispersion of 325 ±25 kms~^ for Gl, with a 
Hernquist profile for the galaxy light distribution. The break radius 
of the galaxy is set to 50 kpc, after which the logarithmic density 
slope becomes 7outer = 3. 

The result from this model is presented in Fig. [14] where the 
inner logarithmic density slope of Gl is shown as a function of the 
external convergence. We see that increasing the external conver- 
gence results in a steeper inner density profile for Gl, as expected. 
It is clear that we need a large mass sheet to explain the differences 
in the slope between using only the VLBI constraints and using the 
joint lensing and dynamics analysis. The most likely value for the 
external convergence from the mass-sheet is /T:ext = 0.65 ±0.15 
(stat.)±0.15 (syst.) to explain this difference, in units of critical 
density. This A^ext for the mass-sheet is much higher than the value 
obtained from just including the group galaxies into the model as 
sub-haloes (A^ext = 0.12 at the position of image A0), which im- 
plies there would have to be a significant underlying parent halo to 
represent the group dark matter distribution. Such a large external 
convergence would be expected to produce further multiple imag- 

^ For the 2SIS+shear+field model. A similar value is obtained at the posi- 
tion of image B. 
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Figure 13. The results of a grid-based modelling analysis of the extended arc emission from the 62108+213 host galaxy. (Left) The observed emission, 
modelled emission, image residuals and source reconstruction for an isothermal mass model. (Right) The same plots for a power-law model with a logarithmic 
density slope of ^' = 2.50. Note that the images are orientated as they were observed with the HST, that is, rotated by 96 degrees east of north. 
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Figure 14. The external convergence of the group halo as a function of 
the logarithmic density slope of Gl. The contours shown are for 68, 95 
and 99 per cent confidence limits. The most likely value of the external 
convergence for a density slope of ^' = 2.45^Q ]^g and a stellar velocity 
dispersion of 325 ± 25 km s-^ is /text = 0.65 ± 0.15 (stat.) ± 0.15 (syst.). 



ing in the field of B2108-I-213. However, none is seen in the HST 
imaging. 



7.5 Possible systematics that could effect the modelling 

The lensing mass models investigated here have been constrained 
by the flux-ratios of the lensed images, leading to the deviation 
from isothermal profiles for Gl and G2. It is worth noting that the 
lensed source is a compact flat- spectrum radio source, which in the 
most luminous cases are typically variable across the whole spec- 
trum. However, the flux-ratios of the two lensed images over optical 
and radio wavele ngths, even when observed at different epochs, are 
almost identical jMcKean et ^1120051 : iMore et al.ll2008h . This sug- 
gests that source variability coupled with a lensing time-delay is not 
causing the flux-ratio between images A and B to differ from what 
is expected from an isothermal model. Furthermore, the lens mod- 
els that use only the extended arc emission from the quasar host 
galaxy, whose emission is not expected to be variable, favour den- 
sity profiles that are steeper than isothermal. The radio spectra of 
the two lensed images between 1.4 and 8.46 GHz are also similar, 
ruling out free-free absorption as a possible cause for the flux ra- 
tio n ot being as expected for a simple isothermal mass distribution 
fe.g. lMittaleTaDl200l l Winn. Rusin & Kochanek|[2003h . 

Another, more likely mechanism that could be possibly re- 
sponsible for cha nging the flux-ratio of the lensed images is mass 
substructure (e.g. IPalal & Kochanekll2002b . An inspection of the 
high resolution //iSr imaging of B2108-I-213 does not show any evi- 
dence for small-scale structure in the form of dwarf galaxies near to 
the lensed images. Such dwarf galaxies have recently been found to 
cause anomalous flux-ratios in some gravitational lens systeni s (e.g. 
iMcKean et al]l200l iMore etaDl2009l : iMacLeod et al.ll2009h . The 
radio emission from the lensed images shows some extended emis- 
sion on VLBI scales, but it is not sufficient to investigate whether 
there is mass substructure near to the lensed images. However, 
the modelling of the lensed emission using the extended gravita- 
tional arc found that the arc structure could be well reproduced by 
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a globally smoo t h mass model with n o significant residuals (e.g. 
lKoopmansl[2005l : IVegetti & Koopmansli2009 a). 

The lensing and dynamics models rely on determining the 
mass enclosed within the Einstein radius, which is dependent on the 
redshifts of the lens and source being known. Although there is no 
doubt over the measured redshift of the lens, the source redshift is 
less secure since it relies on the detection of a break in the spectrum, 
as opposed to several definite detections of emission and absorption 
lines. It is worth noting that if the source redshift is higher than 0.67 
then there will be little change in the enclosed mass due to the cur- 
vature in the angular diameter distances as a function of redshift. If 
the source redshift is lower, then the mass of the lens would need to 
be higher to explain the wide image separation, which would force 
the lensing and dynamics models to produce a flatter logarithmic 
density slope. Therefore, further observations to confirm the source 
redshift should be carried out to remove this possible error in the 
lensing and dynamics models. 



8 DISCUSSION 

B2 108+2 13 is an interesting system because i) the wide separation 
between the two lensed images of 4.56 arcsec suggests a group- 
scale lens and ii) previous attempts to model the system have re- 
quired that the main lensing galaxy has a steeper than isothermal 
density profile. We now discuss to what extent our results have ex- 
plained these unusual features of this lens system and what it tells 
us about the mass distribution of a group-sized halo. 



8.1 Merging groups or a duster? 

We first discuss whether the lens system lies in a group or a clus- 
ter. Our spectroscopic survey of the B2 108-1-2 13 field determined 
redshifts for 90 galaxies. Of these galaxies, up to 52 (including the 
main lensing galaxy Gl) were found to be at redshift ^-^0.365. Thus, 
we have confirmed that there is a large-scale structure associated 
with this gravitational lens system. The total size of the structure is 
between 41 to 52 galaxies, depending on the linking kernel used to 
find the group membership. From Figs.[2land[3]it is clear that there 
is a well defined red sequence in the large-scale structure. This is 
not too surprising considering that we preferentially selected red 
galaxies by their colour. Since we did not search for the less mas- 
sive late-type galaxies during our survey it is likely that the galaxy 
membership of the large-scale structure could be much larger. The 
large number of members already con firmed woul d suggest that 
this structure has an Abell richness of (lAbelil958h . The velocity 
dispersion of the galaxies is between 451 to 694 km s~^ depending 
on the linking kernel used. This is potentially much larger than the 
typical velocity dispersions of low redshift groups, which tend to be 
around 200-450 km s"^ (Zabludoff & Mulchaey 1998), and would 
point toward this structure being at most a poor cluster. However, 
we also found that the velocity distribution of the galaxies has a 
spatial structure, suggesting that there may be two galaxy groups in 
the process of merging to form a cluster. Further observations are 
required to confirm this. 

The X-ray data for this system give the strongest indication 
that this structure is a group of galaxies. The X-ray luminosity of 
the diffuse intra-group gas is Lx,boi,500 1-2 x 10^^ ergs s~^, 
which is toward the low en d of the luminosity distribution for 
clusters at a similar redshift (lEbeling et all |200QI : I Kocevski et al. 
l2007h . Furthermore, the X-ray luminosity is over an order of 



magnitude lower compared to what is expected for intermedi- 
ate redshift groups and clusters with velocity dispersions as large 
694 kms~ ^, the velocity disperion found for the largest link- 
ing kernel C WilHs et al]l2005l : Ijeltema et al.l l200^: 'Fassn acht et all 
2008). Note that for the 450 and 650 kms~^ linking kernels, the 
group velocity dispersion is ^ 470 km s~^, placing the B2 108-1-2 13 
group on the Lx-o't; relation for intermediate redshift groups (see 
lFassnachtet"aDl2008h . The X-ray emission also does not extend 
over the spatial distribution of the galaxies in the structure, but it is 
almost centred on the position of the brightest group galaxy, lens- 
ing galaxy Gl - the projected distance betw een the X-ray peak an d 
lensing galaxy Gl is ^36 kpc (see Fig. [9l Ipassnacht et al.ll2008h . 
Finally, we find that the main lensing galaxy is a cD galaxy, which 
must have formed through the merging of smaller mass structures. 
Therefore, we believe there is strong evidence that this is a dynamic 
multiple-group system in the process of merging. 

There are a number of avenues which can be taken to explore 
the full extent of the B2 108-1-2 13 large-scale structure in the fu- 
ture. The most useful would be to carry out a weak lensing analysis 
from wide field imaging of the system. This would show whether 
the structure has an extended dark matter halo, or possibly shows 
evidence for two or more merging haloes. It would also be interest- 
ing to see whether the dark matter halo follows the X-ray emission, 
or is closer to the emission from th e galaxies - as has been recently 
observed in merging clusters fe.g. lClowe et al.|[200^ : iBradac et al.l 
l2006l.[2008h . Attempts to measure a weak lensing signal from HST 
and ground-based data are currently underway. 

8.2 Mass distribution of a group-scale halo 

Having established that the image- splitting of B2 108-1-2 13 is prob- 
ably caused by a group-scale halo, we now discuss what the gravi- 
tational lensing and stellar dynamics data can tell us about the dis- 
tribution of mass within the structure. 

Including the contribution of the group, either as sub-haloes 
and/or as a mass gradient (or constant mass-sheet) to represent an 
underlying parent group halo, results in lowering the Einstein radii 
of the two lens galaxies between the lensed images, and hence their 
projected mass determined from gravitational lensing. This high- 
lights the importance of carrying out large spectroscopic surveys 
of the environments of lens systems for calculating the projected 
mass of lensing galaxies. The total convergence due to the sub- 
halo population is n ^0.12, which is higher than that found for 
other gravitational lens systems associated with groups and/or hav- 
ing groups found along the line-of-sight (e.g. Fassnacht et al. 2006': 
Momcheva et al. 2006; Auger et al. 2007; Treu et al. 2009). This 
is due to the other gravitational lens systems being part of loose 
groups of only a few members. For B2108-1-213, there are up to 52 
confirmed group members. Also, from Fig. [T] it is clear that there 
is a compact structure of 5 galaxies that are only 15 arcsec from 
the lens system, or at a projected distance of ^75 kpc. Therefore, 
B2 108-1-2 13 probes the mass distribution at the centre of a compact 
structure. 

Our models using only the lensing data obtained from the high 
resolution imaging with VLBI or the HST strongly favour steeper 
than isothermal density profiles at the position of the Einstein ra- 
dius, in the absence of some external convergence field. Such a 
steep density profile for the main lensing galaxy is certainly not 
unique (e.g. Koopmans et al. 2006) and has been predicted by sim- 
ulations for haloes that have recen tly und ergone some kind of inter- 
action ( Dobke, King & FellhaueJEoOTb . However, the steep den- 
sity profile appears to be inconsistent with what was found when 
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the stellar kinematics data are incorporated into the model, which 
requires the average logarithmic density slope between the Einstein 
and effective radius of Gl to be close to isothermal. A possible 
explanation for this contradiction is that the mass distribution of 
the main lensing galaxy can not be simply approximated by a sin- 
gle power-law, but by perhaps, a broken power-law with a knee 
somewhere between the effective radius and the Einstein radius of 
the system. There are two physical mechanisms that could produce 
such a broken power-law, which we now discuss. 

First, the total mass density profile predicted from simula- 
tions is known not to be a single p ower-law at cluster scales (e.g. 
ISand et al.|[2004l : ISand et al.|[20oi : the innermost region is domi- 
nated by the brightest cluster galaxy stars, leading to a total mass 
density with a slope approximately isothermal or steeper. At larger 
radii the inner part of the dark matter halo takes over with a typi- 
cal slope between and 1. At even larger radii the profile becomes 
steeper - in accordance with the predictions of numerical simula- 
tions, reaching again isothermality and then an eve n steeper slope 
of ~ 3 towards the virial radius (e.g. Kneib et al. 2003). The so- 
called break radius, where the dark matter profile is approximately 
isothermal, is of order 100 kpc for clusters. In contrast, at the scale 
of massive galaxies, the smaller break radius and the high concen- 
tration of stars in the centre, conspire to create a total mass density 
profile that is consiste nt with a single pow er law over more than 
two decades in radius (iGavazzi et al .112007b . B2 108-1-2 13 explores 
the intermediate regime of a group size halo, for which we expec t 
the break radius to be of order of tens of kpc|fl ('Maccio et al.ll2007l) . 
Thus, we could interpret the change in the average slope between 
the Einstein and effective radius as a signature of a turnover in the 
mass density profile, qualitatively in agreement with simulations. 
However, to explain the flux-ratio of B2 108-1-2 13 would require the 
turnover to occur close to the Einstein radius, that is, at around 
8 kpc, which seems unlikely for a group-scale halo. More radially 
extended data would be needed to pinpoint the location of a possi- 
ble break, the virial mass of the structure, and the actual slope in 
the outer parts. In a future paper, we plan to add weak lensing to 
our constraints to investigate the full extent of the radial profile. 

A second mechanism that could account for the changing 
slope at the inner part of the halo is galaxy truncation. Simula- 
tions have shown that as galaxies fall into a massive potential, 
like in the cases of groups or clusters, the dark matter is removed 
in the outer parts of the galaxy's halo by tidal s tripping (e.g. 
iBullock et al](200ll : iGao et al.ll2004l : Limousin et a l.i r2009 b) lead- 
ing to a truncated halo. Weak lensing measurements around galax- 
ies in clusters have found evidence for halo truncation, with those 
galaxies closest to th e centre of the clus t er showing the larges t 
levels of truncation (IGavazzi et al.l l2004l : iNataraian et aP l2009h . 
The truncation radius of these galaxies is of order ^ 50-60 kpc 
(iHalkola. Seitz & Pannellal l2007l : INataraian et al.ll2009h . It would 
be reasonable to expect a high level of truncation for the lensing 
galaxy Gl - the galaxy sits close to the centre of the group and is 
likely to have undergone many mergers and interactions over the 
course of time that the group has taken to form. 

A problem with the broken power-law model is that it requires 
the inner part of the halo to be quite flat for the average logarith- 
mic density profile to be close to isothermal. This would essentially 
need Gl to have a more centrally concentrated mass distribution. 
To test this, we considered a model with a constant mass-to-light 
ratio (as opposed to M oc L^^^ from the FJ-relation) for Gl. We 



For our models we assumed a break radius of 50 kpc. 



find that although the lensing constraints can still be fitted with this 
model, the velocity dispersion of Gl is pushed higher, to around 
450 km s~^ and above, since we must now have more mass within 
the effective radius of Gl. This is clearly not consistent with the 
observed stellar velocity dispersion of Gl. A possible solution is 
that there is some external convergence from the group that would 
result in Gl being less massive, but this again would need to be sig- 
nificant (A^ext ^ 0.6) and is not plausible. For this reason, it seems 
that some sort of broken power-law model is not likely. 

Our simplest model for the mass distribution of the lens sys- 
tem that satisfies both the lensing and dynamics constraints, has an 
underlying external convergence gradient. Applying a convergence 
gradient is successful because it changes the expected flux-ratios 
of the lensed images, which from the data currently in hand appear 
not to have been strongly affected by mass substructure or by vari- 
ability. This gives a lens model that is consistent with an isothermal 
mass distribution, and hence with the stellar velocity dispersion of 
Gl. The moderate convergence gradient of only 0.06 arcsec"^ is 
also plausible given the dense environment that the lensing galaxies 
Gl and G2 reside and the gradient seems to increase in the direction 
of the peak in the X-ray emission. Again, weak lensing observa- 
tions of the system would show if there is an external convergence 
gradient across the system, and whether this convergence gradient 
maps out the large-scale structure of the galaxy group associated 
with lensing galaxies Gl and G2. It would also be interesting to 
carry out resolved spectroscopy of the system (e.g. ICzoske et al.l 
l2008h to confirm the redshift of G2 and the lensed quasar, and to 
determine the velocity distribution of Gl to further constrain the 
lens models fe.g. lBarnabe et al.ll2009h . 



9 CONCLUSIONS 

We have presented new optical imaging and spectroscopy of 
the gravitational lens system B2108-I-213. Our observations have 
shown that the main lensing galaxy is the central galaxy of a mas- 
sive group or poor cluster. Using both gravitational lensing and stel- 
lar kinematical analyses, we have carried out extensive modelling 
of the system and have come to the following conclusions. 

i) By targeting a gravitational lens system with a wide image 
separation, we have successfully identified a large group of galax- 
ies at a moderate redshift. Recent studies of the environments of 
gravitational lenses have found many are associated with groups of 
galaxies. However, these systems tend to have lensed image sep- 
arations of order ^ 1-1.5 arcsec, and are therefore, typically as- 
sociated with poor groups, giving only limited information on the 
mass distribution of dark matter haloes in the mass regime between 
galaxies and massive clusters. Further investigations of group-scale 
dark matter haloes should target gravitational lens systems with im- 
age separations similar to or larger than ^ 4 arcsec. 

ii) Our modelling of the inner part of the group mass distribu- 
tion has found that the combined dark and baryonic matter density 
profile is much steeper (i.e. 7' > 2) when compared to what is 
predicted from dark-matter only simulations (i.e. 7 ^ 1-1.5). This 
suggests that in the case of the group associated with B2108-I-213, 
the baryonic component within the inner 20 kpc of the halo centre 
dominates the mass distribution. B2 108-1-2 13 may be a special case 
in this respect since there is a massive cD galaxy in the centre of 
the group. Therefore, it would be useful to find and study group 
systems without a massive central galaxy in the future. 

iii) Our analysis of the gravitational lensing and stellar dy- 
namics data finds that the inner matter density profile is on average 
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isothermal within the Einstein radius, but is steeper (i.e. 7^ ^ 2.5) 
at the position of the Einstein radius. A plausible explanation for 
this contradiction is the presence of an external convergence gradi- 
ent across the lens system. The effect of this gradient is to change 
the flux-ratio of the lensed images and to bring the slope of the 
inner matter density profile determined from lensing to be consis- 
tent with what is predicted from the stellar dynamics of the lensing 
galaxy. The presence of a convergence gradient is certainly plausi- 
ble considering there is a galaxy group associated with Gl and G2. 
Therefore, we find direct evidence for the environment of the lens 
system having an effect on the observed properties of the lensed 
images. Note that by not including a convergence gradient, the in- 
ner matter density profile is artificially steepened and this should be 
considered when modelling other group-scale systems in the future. 

Our conclusions are based on the results from investigating 
only a single group system with gravitational lensing. Clearly, stud- 
ies of many more group systems will be needed to test our find- 
ings for the B2 108-1-2 13 group. Recently, there have been 13 group- 
scale gravitational lens systems identified from the Strong Lensing 
Legacy Survey (S2LS; Limousin et al. 2009a). These groups span 
a large range of group masses and, apparently from maps of the lu- 
minosity distributions, a large scatter of group halo morphologies. 
It will be interesting to see whether this sample of strong lensing 
groups and poor clusters can be modelled without including an ex- 
ternal convergence gradient. 
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